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Intranight Optical Variability of Radio- Quiet Weak 
Emission Line Quasars 
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ABSTRACT 

Based on a recently started programme, we report the first search for intranight 
optical variability among radio-quiet 'weak-line-quasars' (RQWLQs). Eight members 
of this class were observed on 13 nights in the i?-band, such that each source was 
monitored continuously at least once for a minimum duration of about 3.5 hours, 
using the recently installed 130 cm telescope at Devasthal, India. Statistical analysis 
of the differential light curves was carried out using two versions of the J 1 — test. Based 
on the INOV data acquired so far, the radio-quiet WLQ population appears to exhibit 
stronger INOV activity as compared to the general population of radio-quiet quasars 
(RQQs), but similar to the INOV known for radio- loud quasars of non-blazar type. To 
improve upon this early result, as well as extend the comparison to blazars, a factor 
of ^2 improvement in the INOV detection threshold would be needed. Such efforts 
are underway, motivated by the objective to search for the elusive radio-quiet blazars 
using INOV observations. 

Key words: galaxies: active - galaxies: photometry - galaxies: jet - quasars: general 
- (galaxies:) BL Lacertae objects: general - (galaxies:) quasars: emission lines 



1 INTRODUCTION 

Powerful active galactic nuclei (AGN) whose luminosity 
across the electromagnetic spectrum is dominated by a 
Doppler boosted relativistic jet of nonthermal emission are 
termed as blazars. The two subsets of this class, namely BL 
Lac objects (BLOs) and Highly-Polarized-Quasars (HPQs), 
although differentiated by the equivalent widths of emis- 
sion lines, share many properties. But, whereas HPQs have 
an abundant population of weakly polarised quasar coun- 
terparts (mostly radio-quiet quasars, called RQQs), vari- 
ous searches for radio-quiet analogs of BLOs have so far 
remained unsuccessful. BLOs characterized by very weak 
or absent optical/UV emission lines, which have been pur- 
sued in such searche s , are selected from op tical surveys (e.g., 
IJannuzi etHI Il993l : lLondish et~ai1 |2004| ) , alt hough X-ray 
select ed BLOs have also been targeted (e.g., IStocke et al.l 
1990). Usually, the radio-loudness is quantified in terms 
of a parameter R defined as the ratio of the rest-frame 
6 cm to 2500A flux densities and powerful AGN having 
R > 10 are designated as radio-loud (e.g., Jiang et al 



Il989l ). The first radio-quiet AGN showing weak emission 
lines, to be interpreted as a non-BLO was PG 1407+265 at 
z — 0.94, based on the lack of variability on 10 year base- 
line a nd the lack of optical polarization (|McDowell et al.l 
1 19951 : iBerriman et alj rf99Q') . Another example of simi- 
lar spectral peculiarity is the high accretion -rate quasar 
PHL 1811 at z = 0.19 |Leighlv et al.l l2007al ). Samples of 
radio-quiet BLO candidates at low er redshifts (z < 2.2) 
were found in the SPS S survey dYork et al.l |2000| ). by 



ICoUinge et al.l |2005l ) and I Anderson et al.l (|2007l ). and were 
termed 'Weak-Line-Quasars' (WLQs). As a result, dozens of 
WLQs marked by abnormally weak broad emission-lines (i.e, 
rest- frame EW < 15. 4A for the Lv q+NV emission-line corn- 
pie^ Diamond-Stanic et alj 20091) have been reported (e.g., 



120071 : IShen et al.ll201ll ; IStocke et ailll992l : iKellermann et al 
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Although the above studies hav e revealed many WLC 
that are indeed radio-quiet (e.g., IPlotkin et al.l 120101: 
they are commonly identified not as BLOs but RQQs 
having abnormally weak emission lines. This is because, 
in contrast to BLOs (and much like RQQs), the radio- 
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quiet WLQs (RQWLQs) are fo und to exhibit low optical 
polarization (|Smith et al.l 120071 ) and mild optical contin- 
uum variability on t ime scales ranging from days to years 
(|Plotkin et al-feoiObh . This is further supported by the sim- 
ilarity observed between the UV-optical spectral indices, 
a, of W LQs and RQQs. For RQQs t he median value of a 
is -0.52 jDiamond-Stanic et all 120091 ; IPlotkin et alj|2010al). 
as aga inst -1.15 for BLO candidates (e.g., Plotkin et alj 
l2010al ). The reason for the abnormally weak line emission 
in WLQs is yet to be fully understood, but the explana- 
tions proposed basically fall into two categories. One pos- 
sible cause of the abnormality is the high mass of the cen- 
tral BH (Mfljj > 3.1O 9 M ) which can result in an accre- 
tion disk too cold to emit strongly the ioni zing UV pho- 
tons, even when its optical out put is high (|Laor fc Davis! 
l201ll : also, IPlotkin et ail |2010al ). Alternatively, the cover- 
ing factor of the broad-line region (BLR) in WLQs could 
be at least an orde r-of-magnitude smaller com pared to the 
normal QSOs (e.g., iNikoIaiuk fc Waltenl2012i ). An extreme 
version of this scenario is that in WLQs the accretion disk 
is relatively recently established and hence a significant 
BLR is yet to develop (|Hrvniewicz et alJl20ld : lLiu fc Zhaiigl 
l201ll L Conceivably, a poor BLR could also result from the 
weakness of the radiation pressure driven wind when the 
AGN is operating at an exce ptionally low accreti on rate 
(< 10~ 2 to 1 Q- 3 Mj 
lElitzur fc Holl2009h . 

While the above mentioned limited empirical evidences 
and theoretical scenarios are consistent with the quasar in- 
terpretation of the bulk of the WLQ population, they do 
not rule out the possibility of a small subset of the popu- 
lation being, in fact, the long-sought radio-quiet BLOs in 
which optical emission arises predomi nantly from a rela- 
tivist i c jet of synchrotron radiation fe.g.JStocke fc Perrenodl 
Il98ll . iDiamond-Stanic et al.l |2009|, IPlotkin et al.|[2010al and 



Table 1. The 8 RQWLQs studied in the present work. 



e ptionally low accreti c 
picastro et all 120031 
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fe nces therein: also. l stalin fc SrianandB 



One strategy to pursue such a search is to character- 
ize the intra-night optical variability (INOV) of radio-quiet 
WLQs (RQWLQs). It is well established that normal BLOs 
(which are always radio-loud) exhibit a distinctly stronger 
INOV, both in amplitude (ip) and duty cycle (DC), as 
compared t o quasars, specially their radio-quiet majority , 
RQQs (e.g..lGopal-Krishna et al.ll2003l; IStalin et al.|[2004al ; 



iGupta fc Joshill2005l ; ICarini et alil2007l ; iGoval et al.H2012f ). 

From this it is evident that INOV properties can be a strong 
discriminator between blazars and ot her powerful AGN , 
both radio-loud and radio-quiet (e.g, IStalin et al.l l2004al ; 
IGoval et al.ll2012h . The impetus behind our new programme, 
therefore, is to characterize the INOV behavior of RQWLQs 
and the first results are presented here. 



2 THE SAMPLE OF RADIO-QUIET WLQS 

Our sample for INOV monitoring (Table [T]) was derived 
fro m the list of 86 radi o-quiet WLQs published in Table 6 
of (IPlotkin et al.ll2010al') base d on the SDSS Data Release 7 
fDR-7. lAbazaiian et al]|2009t) . Out of that list, we included 
in our sample all 18 objects brighter than R ~ 18.5 which 
are classified as 'high-confidence BL Lac candidate'. Thus 
far, we have been able to carry out intranight monitoring 
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of only 8 of these sources in 13 sessions and the results are 
reported here. 



2.1 Photometric observations 

Continuous monitoring of each RQWLQ was done, mainly 
using the 1.3-m optical telescope (hereafter 1.3-m DFOIq) 
of the Aryabhatta Research Institute of obser vational sci- 
encE S (ARIES), located at Devasthal, India (|Sagar et alj 
l201ll ). DFOT is a fast beam (f/4) optical telescope with a 
pointing accuracy better than 10 arcsec RMS. The telescope 
is equipped with Andor CCD having 2048 x 2048 pixels of 
13.5 micron size, resulting in field of view of 18 arcmin on the 
sky. The CCD is read out with 31 and 1000 kHz speeds, with 
the corresponding system RMS noise of 2.5, 7 e~ and gain 
of 0.7, 2 e~/Analog-to-Digital Unit (ADU). The camera is 
cooled down thermoelectrically to —85 degC. We performed 
continuous monitoring of each source for about 4 hour in the 
SDSS— r passband at which our CCD system has maximum 
sensitivity. For achieving SNR greater than 25-30 our typi- 
cal exposure time was set between 5—8 minutes. The typical 
seeing FWHM during our monitoring sessions was 2 arcsec, 
adequate for these point-like sources. 

One of the RQWLQ (J125219.47+264053.9) was also 
monitored with the 1.04-m Sampurnanand telescope (ST) 
located at ARIES, Nainital (India). Another RQWLQ 
(J090107.60+384659.0) was also monitored using the 2-m 
IUCAA Girawali Observatory (IGO) telescope located near 
Pune (India). The ST has Ritchey-Chretien (RC) optics with 
a f/13 beam (Sagar 1999). The detector was a cryogenically 
cooled 2048 x 2048 chip mounted at the Cassegrain focus. 
This chip has a readout noise of 5.3 e~ /pixel and a gain 
of 10 e~/ADU in the slow readout mode. Each pixel has a 
area of 24 (jm 2 which corresponds to 0.37 arcsec 2 on the sky, 
covering a total field of 13' x 13'. Our observations were car- 
ried out in 2 x 2 binned mode to improve the signal-to-noise 
ratio, and Cousins R filters were used. 

The 2-metre IGO telescope has an RC design with a 
f/10 beam at the Cassegrain focufl The detector was a cryo- 
genically cooled 2110x2048 chip mounted at the Cassegrain 
focus. The pixel area is 15 /1m 2 , so that the image scale of 
0.27 arcsec/ pixel covers an area of 10' x 10' on the sky. 
The readout noise of this CCD is 4.0 e~ /pixel and the gain 



1 Devsthal Fast Optical Telescope 
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is 1.5 e~/ADU. The CCD was used in an unbinned mode 
with Cousins R filters. 

In our sample selection, care was taken to ensure the 
availability of at least two, but usually more, comparison 
stars on the CCD frame that were within about 1 mag of 
the target RQWLQ. This allowed us to identify and discount 
any comparison star which itself varied during a given night 
and hence ensured reliable differential photometry of the 
RQWLQ. 



2.2 Data Reduction 

All pre-processing of the images (bias subtraction, flat- 
fielding and cosmic-ray removal) was carried out using the 
standard tasks available in the data reduction software 
IRAF 0. Instrumental magnitudes of the comparison stars 
and the target source were measured from the frames us- 
ing the Dominion Astronomical Observatory Photometry 
(DAOPHOT II) ifl software designe d for concentric cir - 
cular aperture photometric technique (|Stetsonlll992l . fl987l ). 
As a check on the possible effects of any seeing variations, 
the aperture photometry was carried out with four aperture 
radii, lxFWHM, 2xFWHM, 3xFWHM and 4xFWHM, 
where the seeing disk radius (= FWHM/2) for each CCD 
frame was determined using 5 fairly bright stars on the 
frame. The data reduced using the four aperture radii were 
found to be in generally good agreement. However, the best 
S/N for the differential light curves (DLCs) was nearly al- 
ways found for aperture radii of ~2xFWHM, so we adopted 
that aperture for our final analysis. 

To derive the Differential Light Curves (DLCs) of a 
given target RQWLQ, we selected two steady comparison 
star present within the CCD frames, on the basis of their 
proximity to the target source, both in location and mag- 
nitude. Coordinates of the comparison star pair selected for 
each RQWLQ are given in Table [2] The g — r color dif- 
ference for our 'quasar-star' and 'star-star' pairs is always 
< 1.5, with a media n value of 0.54 (colu mn 7, Table 0. 
Detaile d analyses by ICarini et all (|l992h and Istalin et al.l 
l|2004al ) show that color difference of this magnitude should 
produce negligible effect on the DLCs as the atmospheric 
attenuation changes during a monitoring session. 

Since the selected comparison stars are non-varying, as 
judged from the steadiness of their DLCs, any sharp fluc- 
tuation over a single temporal bin was taken to arise due 
to improper removal of cosmic rays, or some unknown in- 
strumental effect, and such outlier data points (deviating 
by more than 3cr from the mean) were removed from the 
affected DLCs, by applying a mean clip algorithm. In prac- 
tice, such outliers were quite rare and never exceeded two 
data points for any DLC, as displayed in Figure 1. Finally, 
in order to enhance the SNR, without incurring significant 
loss of time resolution, we have taken 3-point box average 
of each DLC. 
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3 ANALYSIS 

Conventionally, the pres ence of INOV in a DLC is quan- 
tified using C-statistics (|jang fc Miller! 119971) . However, re- 
cently Ide Diegol (|2010h has pointed out that this is not a 
valid test as it is based on ratio of two standard devia- 
tions which (unlike variance) are not linear operators and 
the nominal critical value used for confirming the presence 
of variability (i.e., 2.576) is usually too conservative. He has 
therefore advocated more powerful statistical tests, namely, 
the one-way analysis of variance (ANOVA) and the F— test. 
However, a proper use of the ANOVA test requires a rather 
large number of data points in the DLC, so as to have sev- 
eral points within each subgroup used for the analysis; this 
is not feasible for our light curves which typically have only 
around 15 - 20 data points each. Therefore, in this study 
we shall rely on the F— test which is based on the ratio 
of variances as, F = variance( b S erved)/ var i an ce(expected) 
(|de Diegol (|2010l '). Two versions of this test employed in the 
recent literature are: ( i) the standard F— test (hereinafter, 
. (I2012T) and ( ii) scaled F-test, here- 
test Ijoshi et alj (|201ll ). In this work we have 
subjected all our DLCs to both these statistical tests, as 
discussed below. 

An important point to be borne in mind while ap- 
plying the F" 1 — test is that the photometric errors, as re- 
turned by the routines in the IRAF and DAOPHOT soft- 
wares, are normally underestimated by a factor r\ ranging 
between 1.3 and 1.75 , as found in independent studies (e.g. 



F^-test, iGoval et al" 
inafter, F 



Gop al-Krishna et al.ll 19951; ICarcia et al 



2004: 



Stalin et al 



2004bl ; 



5; 
: 



Bachev et al 



.19991 ; ISagar et all 
2005h . In a recent 



analysis of 73 DLCs derived for 73 pairs of 'st eady' stars 
monitored on as many nights, IGoval et al.l (|2012f l estimated 
the best-fit value of r\ to be 1.5. (see, also Sect. 4). The F v — 
statistics can be expressed as: 



pji - 
1 V 2 (° 2 



(q-sl) 



.!>■ 



r) 2 (° 2 s2 } 



J Jr sl~s2 ~ 



u (sl-s2) 



.(1) 



where a 2 q _ 



are the variances 



Dominion Ahtrophysical Observatory Photometry 



and °f s l- S 2) 

of the 'quasar-starl', 'quasar-star2' and 'starl-star2' DLCs 
and {o-g_ s i) = X)£Lo°i.M-r(g - sl)/N, (a 2 _ s2 ) and (a 2 !- a2 ) 
are the mean square (formal) rms errors of the individual 
data points in the 'quasar-starl', 'quasar-star2' and 'starl- 
star2' PLCs, resp ectively. r\ is the scaling factor (= 1.5, cf. 
IGoval et al.ll2012r ). introduced to account for the underesti- 
mation of photometric rms errors returned by the photom- 
etry algorithms used here, as mentioned above. 

The F values computed using Eq. [T] , were then com- 
pared individually with the critical F value, Fc*l lUas , where 
a is the significance level set for the test, and v qs and v sa 
are the degrees of freedom of the 'quasar-star' and 'star-star' 
DLCs, respectively. The smaller the a, the more improbable 
is the result to arise from chance. For the present study, we 
have used two significance levels, a = 0.01 and 0.05, which 
correspond to confidence levels of greater than 99 and 95 per 
cent, respectively. If F is found to exceed the critical value, 
the null hypothesis (i.e., no variability) is discarded to the 
corresponding level of confidence. We have computed sepa- 
rately the F— values for the 'quasar-starl' and 'quasar-star2' 
DLCs (i.e., F? & F$) from Eq.Q] Thus, for a given monitor- 
ing session, a RQWLQ is marked as variable ('V') if for both 
its DLCs F-value ^ F c (0.99), which corresponds to a confi- 
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Table 2. Basic parameters and observing dates of the 8 RQWLQs and their comparison stars. 



IAU Name Date R.A.(J2000) Dec. (.12000) g r g-r 

dd.mm.yy (h m s) (° ' ") (mag) (mag) (mag) 

(1) (2) (3) (4) (5) (6) (7) 



J081250 


79+522531.0 


23.01.2012 


OX 


12 


50 


79 


+52 


2.", 


31 


11 


18 


30 


18 


05 





25 




SI 




08 


12 


50 


27 


+52 


20 


32 


8 


17 


48 


17 


05 





44 




S2 




08 


12 


2!) 


42 


+52 


21) 


4!) 


9 


1!) 


51 


18 


00 


1 


4.3 


J084424 


24+124546.5 


26.02.2012 


08 


44 


24 


24 


+ 12 


45 


40 


5 


18 


2!) 


17 


01 





.37 
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08 


44 


30 


81) 


+ 12 


41 


24 


8 


19 


42 


18 


00 


1 


42 




S2 




08 


44 


39 


26 


+ 12 


44 


54 





18 


27 


17 


8 7 





40 


J090107 


64+384658. 8 


27.02.2012 


09 


01 


(17 


04 


+38 


40 


58 


8 


18 


25 


18 


15 





Of) 




SI 




09 


01 


21 


12 


+38 


42 


14 


1 


18 


93 


17 


00 


1 


24 




S2 




09 


00 


43 


86 


+38 


51 


42 





17 


02 


17 


47 





44 


J090107 


64+384658.8 


16.03.2012 


09 


01 


117 


64 


+38 


40 


58 
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18 


25 


18 


15 





09 




SI 




09 


01 


(111 


15 


+38 


47 


09 


7 


1') 


04 


18 


27 


1 


37 




S2 




09 


00 


59 


94 


+38 


47 


51 


4 


18 


77 


18 


00 





77 
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45+471522.8 
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45 
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15 


22 


8 


17 
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17 
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12 




SI 
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33 
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S2 
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17 
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83 


1 
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45 


+47 
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17 
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17 


53 





12 
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19 


57 


89 
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18 


00 


17 


35 


1 
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S2 
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1!) 
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24 


+47 
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18 
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18 


42 


17 


18 


1 


24 
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47+264053.9 
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12 


52 
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47 
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70 
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52 


37 


93 
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10 
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41 


42 
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10 


71 


10 
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29 
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52 


(111 


81 
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43 


17 


5 


10 


93 


15 


86 


1 


07 


J 142943 


64+385932.2 
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14 


2!) 


43 


64 


+38 


59 


32 


2 


17 


50 


17 


55 





01 




SI 




14 


30 


00 


or, 


+38 


57 


21 


4 


1') 


08 


17 


02 


1 


40 




S2 




14 


29 


30 


69 


+39 


01 


14 


2 


18 


10 


17 


00 


1 


10 


J153044 


08+231013.4 


27.04.2012 


15 


.'ill 


44 


08 


+23 


11) 


13 


4 


17 


83 


17 


59 





24 




SI 




16 


30 


09 


51 


+23 


11 


52 


9 


18 


40 


17 


15 


1 


31 




S2 




15 


30 


47 


70 


+23 


00 


10 


4 


17 


7!) 


17 
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00 
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+23 
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17 
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17 
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15 
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46 
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11 
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1 


17 


2!) 


10 


71 
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S2 




15 


.'ill 


57 
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+2:i 


1)7 


42 


3 


17 


1)1 


16 


05 





30 


J161245 


68+511816.9 


18.05.2012 


16 


12 


45 


68 


+51 


18 


10 


9 


17 


80 


17 


72 





17 




SI 




16 


12 


26 


15 


+51 


22 


14 


6 


18 


1)8 


16 


00 


1 


.39 




S2 




16 


12 


48 


21 


+51 


18 


37 


1 


15 


33 


14 


94 





39 



rms errors on individual points of the 'quasar-star' and 'star- 
star' DLCs, as returned by the DAOPHOT/IRAF routine. 

The threshold criteria for inferring the INOV status of 
a DLC from its computed F-value in this F^-test is identical 
to that adopted above for the f'-test. The inferred INOV 
status of the DLCs of each RQWLQ, relative to two compar- 
ison stars, are presented in Table [3] In the first 5 columns, 
we list the name of the RQWLQ, date of its monitoring, 
telescope used, duration of monitoring and the number, N, 
of data points in the DLCs relative to the two comparison 
stars (si and s2). The next two columns give the computed 
F- values, based on the F v — test and F K — tests. Column 8 
and 9 mention the INOV status of the two DLCs of the 
RQWLQ, as inferred from the F" 7 — test and F K — test, re- 
spectively. Column 10 gives the INOV amplitudes ip derived 
from the two DLCs of the RQWLQ, based on the definition 
given bv lRomero, Cellone. fe Combil (1 19991 ): 

ij) = \/(D max - D min ) 2 - 2o 2 (4) 

with D mintmax — minimum (maxim um) in the RQWLQ 
DLC and a 2 = rfi^l-s), where, 77 =1.5 (|Goval et al.ll2012l ). 
Column 11 lists the square root of the scaling factor, k (Eq. 
3). which has been used to scale the variance of the star— star 
DLCs while computing the F-value in the scaled F-test (Eq. 
2). The last column gives our averaged photometric error 



dence level 99 per cent; non-variable ('NV') if even one of 
the two DLCs is found to have F-value less than F c (0.95). 
The remaining cases are termed as probably variable ('PV'). 

An alternative approach to quantify t he IN OV status of 
a DLC has been followed in I Joshi et al.l l|201ll ). the 'scaled 
F— test'. Instead of rj, this test relies on a factor k equal to 
the ratio of the mean square rms errors of the data points 
in the quasar DLC relative to a comparison star and in the 
DLC of that star relative to the other comparison star. This 
parameter is intended to correct for any bias which may 
arise due to some systematic difference between the pho- 
tometric errors of the data points in the 'quasar-star' and 
'star-star' DLCs (e.g., due to a brightness mismatch between 
the quasar and the comparison star(s)). Thus, in this 'scaled' 
F-test, 



Ff = 



var(q — si) 
k x var(sl — s2) 

with k, defined as, 



Ff = 



var(q — s2) 
k x var(sl — s2) 



(2) 



EN 
i=o CT 



2 

i,err 



(q - a) /N 



£™ < err ( S l-s2)/7V_ 
where Oi. err {q — s) and o~i errisl 



(3) 



<^l- s 2> 

s2) are, respectively, the 
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o"i,err(9 — s) in the 'quasar— star' DLCs (i.e., mean value for 
q-sl and q-s2 DLCs), which typically lies between 0.01 and 
0.02 mag. 



3.1 The INOV duty cycle (DC) 

To recapitulate, a RQWLQ in a given session is marked as 
variable ('V') if its DLCs relative to the two comparison 
stars, are both found to have f-value ^ F c (0.99), which 
corresponds to a confidence level ^ 99 per cent; non-variable 
('NV') if even one of the two DLCs is found to have F- 
value less than F c (0.95). The remaining cases are marked as 
probably variable ('PV'). 

The duty cycle of INOV was computed using the defi- 
nition by I Romero. Cellone. fc Combil (| 19991 ). 



DC 



100 rovM) percent 



(5) 



where Ati = Ati i0 f, s (l + z)^ 1 is duration of the monitoring 
session of a source on the i th night, corrected for its cos- 
mological redshift, z. Since the duration of the observing 
session for a given source differs from night to night, the 
computation of DC has been weighted by the actual moni- 
toring duration Ati on the i th night. Ni was set equal to 1, 
if INOV was detected (i.e., 'V'), otherwise Ni was taken as 
zero. 



4 DISCUSSION AND CONCLUSIONS 

The present study marks the beginning of a systematic in- 
vestigation of the INOV properties of radio-quiet weak-line- 
quasars (RQWLQs). This initial attempt is based on a mod- 
est size sample containing 8 RQWLQs, for which the derived 
results are presented in Table 3. Using the F K — test we ob- 
tained an INOV duty cycle (DC) of ~ 13 per cent which rises 
to ~ 30 per cent if the two cases of probable INOV ('PV') 
are included. On the other hand, the F v — test yields for the 
same dataset an INOV DC of ~ 6 per cent (taking the best- 
fit value of rj = 1.5, Sect. 3). Thus, taken together, the two 
F— tests lead to an average INOV DC of around 9 per cent 
for RQWLQs, for monitoring sessions lasting > 3.5 hours. 
In order to assess the effect of possible uncertainty in the n 
factor (Sect. 3), we have repeated the F v — test for the entire 
sample, taking two extreme valu es of rj (= 1.3 and 1.75), as 
reported in the literature (e.g.. iGopal-Krishna et al. ll 19951 ; 
iGarcia et ailll999l ; IStalin et al1l2004al ; iBachev et al.ll2005l ). 
The computed INOV duty cycles for both these extreme 
values of n are still 6 per cent, i.e., the same as that es- 
ti mated above taking r/ = 1.5, the best-fit estimate given 
in iGoval et""aH (|2012n . Thus, the i^-test is found to give 
consistent results over the maximum plausible range of un- 
certainty in 77. 

At this point it seems worthwhile to also mention the 
DC estimate based on the more conservative, but hitherto 
much more extensively used C— test (Sect. 3). We find that 
the only change to Table 3, resulting from the applica- 
tion of C— test to our dataset is that INOV status of the 
WLQ J121929.45+471522.8 on 26.02.2012 changes from 'V 
to 'PV. This leaves no clear incidence of INOV detection in 
the present data. Treating 'PV cases as 'V yields an INOV 
duty cycle of ~ 6 percent, which would clearly be an upper 



limit, albeit using a small sample. Our subsequent discussion 
will only be based on the results obtained from the F— test 
as it is believed to be a more powerful test (Diego 2010; Sect. 
3). 

Bearing in mind the modest size of our RQWLQ sam- 
ple at this stage, we now attempt a comparison of the INOV 
duty cycle with the estimates available for RQQs and other 
AGN classes, such as non-blazar type flat-spectrum radio 
quasars (FSRQs) and blazars. INOV duty cycles for these 
AGN classes have been extensively reported in the literature 
(e.g., IStalin et all l2004bl ; IGoval et all l2012h . mostly based 
on DLCs longer than ~ 4-hours (which broadly holds even 
for the present DLCs of RQWLQs, as well). One limita- 
tion encountered in making the comparison is that for the 
observations of all these other AGN types, an INOV de- 
tection threshold (ipum) of 1% to 2% had typically been 
achieved (at least in our programme from ARIES, Sect. 1). 
Being 1 to 2 magnitudes fainter, the INOV detection thresh- 
old reached for the present sample of RQWLQs is less deep 
(tpiim ~ 4 — 5%, Table 3). Thus, for the purpose of compari- 
son with the afore-mentioned other AGN types, our present 
estimate of INOV DC for RQWLQs (~ 9 per cent) must be 
treated as a lower limit. It would be very interesting to check 
if a factor of 2 — 3 improvement in ipu m would lead to a much 
higher INOV DC for RQWLQs, perhaps even approaching 
the level of ~ 50% which is established for strong INOV 
(i.e., tp > 3%) of blazars (BL Lacs and high-polarization ra- 
dio quasars) when they are monitored for > 4 hours, (e.g., 



Gopal- Krishna et al.ll2003l; | Stalin et al.ll2004c]|bl;ISagar et all 
20041 ; IGopal-Krishna et al lhoilMGoval et aU2012r ). The DC 
for strong INOV is found to be only ~ 7% for non-blazar 
type FSRQs (based on the f-test, (e.g jGoval et al.ll2012h 
and practically zero for radio-quiet quasars si nce they are 
not known to show INOV with tjj > 3% (e.g., Goval et al.l 
120071 ; IStalin et al.ll2004al . 120051 ; IGopal-Krishna et al.lbOoi T 
Thus, one indication emerging from this first INOV obser- 
vations of radio-quiet WLQs is that their INOV level, as 
a class, is likely to be significantly stronger in compari- 
son to the general population of radio-quiet quasars and, 
indeed similar that that known for non-blazar type radio 
quasars (FSRQs). It remains to be seen whether on attain- 
ing a matching INOV detection threshold ipu m ~ 1 — 2%, 
the INOV activity level of RQWLQs will be found to be 
stronger, perh aps approaching th e high levels exhibited by 
blazars ((e.g., I Goval et al.l |2012| ) and references therein). 
This remains an outstanding question to be pursued, in 
view of its potential for unravelling the nature of WLQs 
and for the key question whether radio-quiet BL Lacs at all 
exist (Sect. 1). It may be noted here that a hint that, com- 
pared to normal RQQs, RQWLQs may show stronger opti- 
ca l/UV variability on year-l ike time scale, has been reported 
bv lStalin fc Srianandl j2005l ); though it is based on monitor- 



ing of just one RQWLQ (SDSS J153259. 96-003944.1 at z 
= 4.67). 

To summarize, the twin objectives pursued in this ex- 
ploratory, first INOV study of radio-quiet weak-line-quasars 
(RQWLQs) are: (a) To find cases of strong INOV (ip > 3%), 
any such RQWLQs would be outstanding candidates for the 
putative radio-quiet BL Lacs, and (b) To quantify the INOV 
duty cycle for the class of RQWLQs, for both strong and 
weaker INOV. In our program we have so far been able to 
cover only a modest size sample containing 8 RQWLQs, each 
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monitored in at least one session lasting > 3.5 hours. This 
has led to the result that the duty cycle of strong INOV in 
this class of AGNs seems to be higher than that known for 
radio-quiet quasars and is similar to that known for (non- 
blazar) FSRQs. This early indication provides impetus to 
continue this programme, in particular, to check if blazar- 
like INOV levels occur in some RQWLQs. To attain the 
required observational capability, a factor of > 2 improve- 
ment in the INOV detection threshold would be needed and 
we are attempting to achieve this by monitoring relatively 
bright RQWLQs on dark nights, possibly using a telescope 
larger than the newly installed 1.3-metre DFOT used in the 
present work. 
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Figure 1. Differential light curves (DLCs), after three point box average, for the 8 RQWLQs in our sample. The name of the quasar 
along with the date and duration of the monitoring session are given at the top of each panel. In each panel the upper DLC is derived 
using the two comparison stars, while the lower two DLCs are the 'quasar-star' DLCs, as defined in the labels on the right side. Any likely 
outlier point (at > 3cr) in the DLCs are marked with crosses (see Sect. 2) and those points are excluded from the statistical analysis. 
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Table 3. Observational details and INOV results for the sample of 8 RQWLQs. 



RQWLQ Date Tel. T F-test values INOV status^ INOV amplitude yfi£_ yj {a-^ err ) 

dd.mm.yyyy used hr F\ ,F$ F£,F£ F^-test F K -test V>l(%)> (q-s) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 
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.57 





.02 



a The number of data points after three point box average. 

^ V=variable, i.e., confidence level ^ 0.99; PV=probable variable, i.e., 0.95 — 0.99 confidence level; NV=non-variable, i.e., confidence level < 0.95. 

Variability status values based on quasar-starl and quasar-star2 pairs are separated by a comma. 

c Here k = ("'q_ s )/(^J 1 _ s 2) ( as m Eq.|3jl, is used to scale the variance of starl-star2 DLCs for the scaled F-test. 



